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1.0 GENERAL INFORMATION 


1.1 As a result of the combined efforts of the U.S. Naval Oceanographic 
Office (NAVOCEANO) and the British Geological Survey (BGS), the new Department 
of Defense (DOD) magnetic model, WC-85, is available. WC-85 has been adopted 
as the official model by United States and United Kingdom defense 
establishments, and by the International Hydrographic Bureau. The U.S. models 
the Earth's main magnetic field; the U.K. models the secular variation or the 
slow changes in the main field over time. John M. Quinn of NAVOCEANO and 
David R. Barraclough and David J. Kerridge of BGS coordinated the WC-85 
modeling effort. 


1.2 The WC-85 model is defined by a set of spherical harmonic coefficients to 
degree and order 12 which describes the main magnetic field of the Earth at 
the base epoch of 1985.0 and by a set of spherical harmonic coefficients to 
degree and order 8 which predicts the slow annual changes of the Earth's main 
magnetic field for 1985.0 to 1990.0. The main field coefficients at any epoch 
other than the base epoch are adjusted linearly with time using the secular 
variation coefficients. The model gives a global description of the northward 
(X), eastward (Y), and downward vertical (Z) magnetic field components, total 
and horizontal intensities (F and H), the declination (D) and the inclination 
(I) of the geomagnetic field for altitudes to 800 km. Because the Earth's 
magnetic field changes with time, the model is updated every five years. 
Thus, WC-85 supersedes WC-80 and will be superseded by WC-90. 


1.3 The WC-85 main field model is based on magnetic-vector observations from 
the MAGSAT satellite and from NAVOCEANO's Project MAGNET aircraft surveys. 
The secular-variation model is based on historical data and annual means of 
ae al observations from worldwide magnetic observatories (Quinn 
et al., 1986). 


1.4 The WC-85 model represents a very smooth interpretation of the Earth's 
main magnetic field and can resolve only core-related features with 
wavelengths of about 3500 km or longer. Therefore, model values may not be 
sufficient for some uses in surface and air navigation where declination 
anomalies and local magnetic features appear (Vega and Jack, 1986). 


2.0 EXPLANATION 


2.1 The World Magnetic Model describes’ the eomagnetic potential 
V(r, 6, 6, t) by the following mathematical expression (ouine. 1986) : 


12 n 
V(r, @, 6, t) =a »; »: (r/a)"*? {a™(t) cos mi + hn(t) sin mh pn (cos 6), 


n=1 m=0 
where in spherical coordinates a is the mean radius of the Earth (6371.2 km), 


r is the radial distance from the Earth's center, @ is the colatitude, ¢ is 
the longitude, t is the time, ph(cos 6) are the Schmidt quasi-normalized 


, 


LSA I SOS 


associated Legendre functions, and gq, and hr are the spherical harmonic : 


coefficients. The secular variation of the geomagnetic potential is the time 
derivative of the above expression truncated to degree and order 8. The 
northward (X), eastward (v) and vertically down (Z) components of the magnetic 


induction Blr r, 0, +, t)= -VV are given as follows (Quinn et al., 1986): 


X = -B, = (1/r) aV/a0; 
Y = By = -(r sin a)” -1 OV/d; } 
Z = -B. = ev/or e " 


2.2 The ASCII FORTRAN Subroutine GEOMAG is available that computes the 
elements of inclination, declination, and total intensity for any geodetic 
position above the International Astronomical Union (1966) reference ellipsoid 


and epoch within the 5-year lifespan of the model. In geodetic coordinates D, " 

I, F, and H are related to the X, Y, and Z components of the magnetic u 
induction as follows: 

Y 

pD = tan! (yx) Fe Vx? + vy? + 7° ‘ 

1 = sin7! (z/F) H= Vx2 4 y2 , ’ 


The program uses any set of Schmidt normalized spherical harmonic « 
coefficients, not to exceed degree and order 12, such as IGRF-80, DGRF-75, : 
WC-85, and others. 


2.3 Tables of grid values of X, Y, Z, D, I, H and F at every 5-degree grid “* 
intersection of latitude and longitude for 1985.0 at the Earth's surface are 
provided. Annual change values (i.e., main field values at the base epoch 
plus one year minus main field values at the base epoch) are given as the 
second entry for each point. The D and I values are given in degrees decimal 
and their rates of annual change in minutes per year. Positive values equal 
east declination and downward inclination; negative values equal west ‘t 
declination and upward inclination. X, Y, Z, H, and F are in units of ul 
nanoteslas for the main field values and nanoteslas per year for the secular 
variation values. Positive values for Z are taken downward; negative values 
are taken upward. For example, WC-85 model values at 10°N and 140°E are as ‘ 


follows: 
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f: Main Field Annual Change 1 

? ( 

x 37039 nT 4.5 nT yr=} v 

Y 1092 nT - 2.5 nT yro- ! 

Z 2308 nT -14.3 nT yyol : 
D 1.7° bad e2! yr" 

I 3.6° - 1.3! yest x 

H 37055 nT 147 yr~ ¥ 

: F 37127 nt 25 nT yrol ‘ 

iE : 

| Small-scale magnetic contour charts of X, Y, Z, D, I, H, and F, and of their ‘ 

f corresponding annual changes are also provided. The map projection is 4 

te Mercator. ; 

Pa! 

ie 2.4 The main field model for 1985 gives the following positions for the "4 

t magnetic dip-poles at 1985.0 (Quinn et al., 1986): ‘: 

F North dip-pole 77.5°N 102.7°W | 


South dip-pole 65.2°S 139.2°E 


TPP ae al 


2.5 For more information on the availability of the WC-85 coefficients and 
GEOMAG subroutine on magnetic tape, write to: 


Commanding Officer n 
U.S. Naval Oceanographic Office ; 
Bay St. Louis 

NSTL, MS 39522-5001 


Attention: DOD Geomagnetic Data Library 4 
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TABLE 1. SPHERICAL HARMONIC COEFFICIENTS - WC-85. 


MAIN FIELD SECULAR CHANGE 


1 0 -29879.8 .0 21.9 .0 
Ye 1 1 -1903.3 5490.5 10.6 -31.5 
2 0 -2070.6 .0 ate .0 
- 2 1 3040.8 -2189.1 1.8 -9.7 
te 2 2 1696.7 -312.0 9.3 -19.9 
3 0 1303.9 .0 8.3 .0 
3 1 -2203.0 -310.3 27.6 6.1 
fe 3 2 1241.7 282.6 -.6 1.3 
s 3 3 839.4 -299,2 2.4 -13.0 
4 0 933.8 .0 “fe2 .0 
2 4 1 781.8 227.2 wl 1.3 
az 4 2 359.0 -246.7 -9.7 3.6 | 
4 3 -424.5 72.5 -1.7 2.5 
on 4 4 164.5 -299.1 -9.3 .6 
f 5 0 -216.4 .0 1.4 .0 
5 1 353.0 43.4 -.5 -.9 
r 5 2 254.3 148.2 -1.2 .6 
7 5 3 -93.7 -154.8 =22 3 
5 4 -157.5 -71.8 9 2.4 
5 5 -45.2 91.5 .0 -1.4 
Gs, 6 0 53.2 .0 3.1 .0 
ne 6 1 63.8 -12.3 .0 af 
6 2 51.3 87.9 1.8 2.1 
6 3 -188.4 67.8 52 -1.4 
) 6 4 3.3 “51.1 -.4 -4.3 
4 6 5 20.3 -4.0 2.4 a 
6 6 -101.7 20.8 1.8 .0 
FE 7 0 76.9 .0 ats | .0 
= 7 1 -60.7 -80.1 -.8 .0 
7 2 Py -25.9 #1.2 1.2 
7 3 25.4 -.9 1.1 2.0 
FE 7 4 -8.1 21.6 0 206 
7 5 6.9 18.5 .6 9 
7 6 7.0 -20.0 -1.8 8 
u 7 7 “4.4 “7.7 -1.2 4 
" 8 0 18.4 0 22 0 
8 1 5.1 3.8 .0 -.6 
% 8 2 1.2 -20.2 a -1.5 
ng 8 3 -12.0 5.0 ol el 
8 4 -9,1 -24.2 02 -1.1 
, 8 5 ol 12.2 -.3 04 
: 8 6 4.7 7.6 -el -2.0 
8 7 6.5 -16.3 22 9 
. 8 8. -9.5 -10.9 -2.2 1.5 
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SPHERICAL HARMONIC COEFFICIENTS - WC-85 


MAIN FIELD 


° ° e e 


e 
NORM OW OK UN HAUP ANP NOWAAUOHY POO 


° 
Ol ell A od 
eee e @ @ 


7 
9 
9 
2 
5 
3 
0 
5 
5 
8 
4 
7 
5 
5 2 
1 5 
6 4 
2 a 
6 -1 
9 3 
8 is 
2 6 
3 
8 
0 
1 
2 
4 
4 
6 
5 
7 
3 
5 
8 


mW tt eet ) 


oe e 
Wm © 


(Con.). 


SECULAR CHANGE 


Qqooocncnocoeo°o°ee°oc & 


Sim sss ter 


¥is 2, ie | PPUPLS | 


coe = fd po | 


ae ne |) oe a 


TABLE 2. WC-85 MAIN FIELD AND ANNUAL CHANGE GRID VALUES 
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